Crossover from time-correlated single-electron tunneling to that of Cooper pairs 



o : 
o . 

(N ■ 

> 
O ■ 

^: 
i> : 

"a' 

o ' 
o 

u 

^: 
^ : 

^ ■ 

I 

o : 
> 

(N 

o : 
o . 



X 



Jonas BylanderEI Tim DutyQ Goran Johansson, and Per Delsing 
Microtechnology and Nanoscience, MC2, Chalmers University of Technology, SE-412 96 Goteborg, Sweden 

(Dated: February 1, 2008) 

We have studied charge transport in a one-dimensional chain of small Josephson junctions using a 
single-electron transistor. We observe a crossover from time-correlated tunneling of single electrons 
to that of Cooper pairs as a function of both magnetic field and current. At relatively high magnetic 
field, single-electron transport dominates and the tunneling frequency is given by / = //e, where 
I is the current through the chain and e is the electron's charge. As the magnetic field is lowered, 
the frequency gradually shifts to / = 7/2e for I > 200 fA, indicating Cooper-pair transport. For the 
parameters of the measured sample, we expect the Cooper-pair transport to be incoherent. 

PACS numbers: 73.23.Hk, 73.23.-b, 85.35.Gv, 85.25.Cp 



Charge transport in ID and 2D arrays of small 
Josephson junctions exhibits a wide range of physical 
phenomena.— In these systems there is a competition 
between the Coulomb blockade, which tends to localize 
charge, and the Josephson effect, which tends to delocal- 
ize it. Depending on the parameters of the Josephson 
junctions in the array, the transport can be described in 
terms of either vortices or charges, which are dual enti- 
ties in a superconducting system in the sense that phase 
and charge are conjugate variables. For strong Josephson 
coupling, Ej , the transport is better described in terms of 
vortices. On the other hand, if the charging energy, Ec, 
is larger, the system is better described in terms of charge 
transport. This duality is not perfect since charge can be 
carried by either Cooper pairs or electrons, whereas there 
is only one type of vortex. Thus the competition between 
Cooper-pair tunneling and single-electron tunneling is of 
particular interest. 

We have previously demonstrated that single electrons 
can be counted one by one as they tunnel through a 
ID series-array of small metallic islands connected by 
Josephson junctions.^ As one excess electron charges an 
island in the array, it polarizes the neighboring islands 
and forms a single-charge "soliton."'^-'' Different solitons 
affect each other by Coulomb repulsion, and therefore 
they form a ID Wigner-like lattice that moves along 
the array. This spatial separation enables a detector to 
resolve the individual charges as they pass by. More- 
over, their passage is time correlated^j^ii with the fre- 
quency / — I/e, where / is the current and e the elec- 
tron's charge. As detector, we used a single-electron 
transistor*'^ (SET) connected to the end of the array. 

In Ref. '3 we only discussed the single electron trans- 
port. However, since the array is superconducting, the 
current can be carried either by electrons or by Cooper 
pairs. In this paper, we report new results from mea- 
surements on the same device (Fig. [Ij, where we now 
study the competition between single-electron tunneling 
and single Cooper-pair tunneling. We demonstrate a 
crossover from single-electron transport to Cooper-pair 
transport as a function of magnetic field and current. 

The studied array is in the strong charging limit, 
Ec S> Ej , and is consequently best described in terms of 



charge transport. Apart from single electron tunneling, 
we can in principle have two different kinds of Cooper- 
pair tunneling: coherent and incoherent. The former 
is equivalent to Bloch oscillationSfiSiii where the sys- 
tem adiabatically follows the lower energy band of each 
junction without dissipation. The latter involves transi- 
tions to excited states and exchange of energy with the 
environment li^iii^ In this particular sample, Ej is smaller 
than the thermal energy ksT, and the system can there- 
fore easily be excited, leading to dissipation and loss of 
coherence. Thus any Cooper-pair transport in our device 
should be predominantly incoherent. 

We fabricated the sample using e-beam lithography 
and triple angle evaporation of aluminum, allowing us 
to use different oxidation parameters for the array and 
SET junctions. The average normal state resistance for 
each of the = 50 array junctions was 940 kO, which 
gives Ej/ks ~ 10 mK. The charging energy per junc- 
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FIG. f : (a) Scanning electron micrograph of the sample. The 
last few islands of the 50-junction array are shown. The SET's 
source, drain and gate electrodes are labeled S, D and G, re- 
spectively, (b) Simplified circuit diagram. The charge enter- 
ing the SET island modifies the dissipation of the LC res- 
onator, which is detected by RF reflectometry. 
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tion was Ec/kB = 2.2 K, corresponding to a junction 
capacitance of C = 0.42 fF. The capacitance to ground 
of each island was Cp » 0.03 fF, giving a single-electron 
soliton size A « \/cJCq ~ 4, which is the number of is- 
lands over which the array is polarized by a single excess 
charge. The SET source-drain resistance was SOkfJ. A 
SEM picture of a sample is shown in Fig.[TJa), and ar- 
ray current-voltage characteristics in the lower graph of 
Fig.H 

We performed the measurements in a dilution refriger- 
ator at approximately 30 mK. A magnetic field of up to 
3 T could be applied parallel to the substrate. We deter- 
mined the parallel critical field for our sample from I —V 
curves of the SET to be ^n^^ « 650 mT. 

The SET was embedded in an LC circuit and operated 
in the radio- frequency mode (RF-SET);— the resonance 
frequency was 358 MHz and the bandwidth 10 MHz. The 
circuit's reflection coefficient depends sensitively on the 
charge induced on the SET island. After amplification 
by cold and room temperature amplifiers, the reflected 
signal was demodulated by homodyne mixing and the 
baseband signal was then measured by a spectrum ana- 
lyzer, see Fig.[TJb). The charge sensitivity is, in general, 
magnetic field dependent, but was typically 20 /ie/\/Hz 
in our measurement. 

The array was biased using a Keithley 263 Calibra- 
tor/Source in feedback mode to maintain a constant av- 
erage current. The biasing line for the array was heavily 
filtered using both stainless steel powder filters and com- 
mercial filters.— 

When a constant bias is applied, charge solitons move 
through the array and approach the SET. The space cor- 
relation of the Wigner lattice translates into time cor- 
related tunneling of charges into the SET at the end of 
the array. Since the full tunneling charge is injected into 
the SET island, the SET acts as a non-linear charge de- 
tector. Numerical simulations'^ show, that the charging 
of the SET island occurs quasi-continuously, whereas the 
discharge happens abruptly by a tunneling event. Be- 
cause of the limited bandwidth of our detector, we can 
only follow the gradual charging, but not the much faster 
tunneling event. 

Our limited sensitivity prevents us from discriminating 
the gradual charging due to single electrons from that of 
Cooper pairs, as either tunneling event gives rise only 
to a "click" in the detector response. However, we can 
discriminate by frequency; if the current is carried by 
electrons, the frequency will be fe — I /e, whereas if it is 
carried by Cooper pairs it will be /2e = //2e. The power 
spectrum of the signal will thus reveal information about 
the type of charge carrier. 

The upper graph in Fig. [5] shows the power spectral 
density of the output signal from the mixer for several 
different values of magnetic field, when the array is biased 
with a constant current of 275 f A (from top to bottom, B| | 
goes from 100 to 500 mT in steps of 50 mT). For a field of 
500 mT the spectrum has a clear peak at the frequency fe 
(the dashed line to the right in Fig. [21). This peak is due 
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FIG. 2: (color online) Upper graph: Power spectral den- 
sity of the output signal from the RF-SET when the cur- 
rent through the array is maintained at 275 fA. The curves 
have been displaced vertically for clarity. Black lines are 
data starting at i3| | = 500 mT (bottom) and continuing ev- 
ery 50 mT to B|| = 100 mT (top). The solid (red) lines 
are fits to a Lorentzian plus a ~ 1// background. The two 
dashed (blue) lines correspond to / = J/e = 1.72 MHz and 
/ — I /2e = 0.86 MHz. Lower graph: Array I — V curves at 
B|| = (black), 400 (red), and 800 mT (blue). 



to time correlated transport of single electrons.^ At fields 
higher than 500 mT, the I — V characteristic becomes 
very steep, and thus the current becomes very sensitive 
to fluctuations in the bias and to background charges, as 
discussed in Ref. |2|. Therefore, the peak in the spectral 
density is smeared and disappears into the noise floor. 
For decreasing magnetic field, the peak gradually moves 
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to lower frequencies, and around 200 mT it appears at f2e 
(left dashed line in Fig. (2]). At even lower fields, the peak 
is smeared and could not be observed below 100 mT. 

We define an effective charge as the nominal array cur- 
rent divided by the peak frequency, QeS — I / fpcak, as 
obtained from fitting to a Lorentzian and a 1//" back- 
ground, where a ^ 1. In the intermediate regime, where 
le < Qoff < 2e, there is a mixture of extra single elec- 
trons and Cooper pairs in the array. In Fig. EJa), we 
show how QcS changes as a function of magnetic field 
for a fixed bias current / = 200 fA. For < 250 mT, 
Cooper pair transport dominates; in the intermediate 
regime, 250 mT < < 400 mT, there is coexisting le 
and 2e transport; and for > 400 mT, there is pre- 
dominantly single-electron transport. In Fig. [3l[b), we 
show how Qcff varies as a function of current for a fixed 
magnetic field = 150 mT. Here, the mixed le and 2e 
transport occurs in the region below / « 200 fA, whereas 
Qcff = 2e above this current. 

Figure IHa) shows that this magnetic field-induced 
crossover occurs only for relatively high current; at low 
current, Qoff — le for all magnetic fields. Moreover, 
this figure shows that the current-induced crossover oc- 
curs only at low magnetic field, where low current favors 
electron transport whereas high current favors Cooper- 
pairs. The measured voltages (in mV) across the array 
are shown as contours. 

In Fig. HKb), we display the normalized width of the 
peak (half width divided by frequency) for the same cur- 
rents and magnetic fields as in (a). It is clear that the 
sharpness, i.e., the degree of correlation between suc- 
cessive tunneling events, is greater when the transport 
is dominated by single-electron tunneling, and for small 
currents where there are few solitons inside the array at 
a given time. 

The fact that charge transport with only Cooper pairs 
is less correlated than transport with quasiparticles can 
be qualitatively explained using energy arguments. The 
real part of the impedance seen from a junction inside 
the array is much smaller than the quantum resistance 
Rq = /i/4e^ ~ 6 kfi ^ why energy exchange with the envi- 
ronment is very ineffective for Cooper pairs. Therefore, 
they tunnel in principle only when the charging energy 
difference before and after the tunneling event is smaller 
than Ej. On the contrary, quasiparticles are energeti- 
cally allowed to tunnel as soon as the charging energy 
difference is positive. Thus, the Cooper pairs are more 
prone to be trapped inside the array, degrading the time 
correlation. Fluctuating background charges, and also 
approaching solitons, can change the local bias of a junc- 
tion so that the elastic channel opens and the Cooper pair 
tunnels. In the regime with mixed le and 2e transport, 
the peak sharpening when Qcff approaches le suggests, 
that the more mobile le solitons (quasiparticles) are in- 
deed effective in "freeing" the 2e solitons (Cooper pairs). 

The situation is thus rather complex with a number of 
things that affect the transport, including randomly dis- 
tributed background charges, non-equilibrium quasipar- 
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FIG. 3: (color online) (a) Effective average charge Qotr = 
///peak vs. parallel magnetic field _B|| for the current J=200fA 
through the array. For high field, the charge is le, whereas at 
lower field it increases and reaches 2e around 200 mT. Inset: 
Magnetic field dependence of the threshold voltages for single 
electrons (Vf,e) and Cooper pairs (Vt,2e), see Eq. ((T)). In the 
regions labeled "le" and "2e" only single electrons and single 
Cooper pairs are allowed, respectively. Above both thresh- 
olds, both types of charge carriers are allowed, (b) EH'ective 
charge versus current for = 150 mT; see discussion in the 
main text. The peak frequencies /peak and the error bars 
in both plots are obtained from fitting each power spectrum 
to a Lorentzian, see Fig. 2, and Qeff is calculated using the 
nominal array current I. 



tides of unknown density, and the electromagnetic envi- 
ronment. Therefore, a complete quantitative description 
of these results is hard to attain. We can, however, give a 
number of qualitative arguments to explain the observed 
phenomena. 

Which type of transport will dominate is largely de- 
termined by the type of carrier that is being injected at 
the first junction, since well inside the array, the charges 
repel each other. The threshold Vt^e for injecting a single 
electron depends on the magnetic field since the energy 
2A has to be supplied to break a Cooper pair. At zero 
temperature, Vt can be calculated from electrostatic en- 
ergy considerations.^^ In our limit 1 <C A^ <C N'^ (and 
ignoring the effect of random background charges), we 
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FIG. 4: (color online) (a) Effective charge versus magnetic 
field and current. The dashed (blue) lines indicate the cut- 
outs shown in Fig. 3 for fixed field and current, (b) Ifalf width 
of the peak scaled by peak frequency. The solid (green) lines 
and labels in both plots are contours of constant bias voltage 
in mV across the array. 

get for electrons and Cooper pairs, respectively: 
y^B) = ^[l + exp(-l/A)] (l + 

^*.2e=22^[l+CXp(-l/A)]. 

Here CeS = a/C'o + 4CCo = 0.23 fF is the effective island 
capacitance, Ec — e^/{2C + Cq + Ccs) — l.TKfcs is the 
first island's charging energy, and A(0)/fcs = 2.4 K is the 
superconducting energy gap of our aluminum thin films 
at zero magnetic field and temperature. The dependence 
(H]) is displayed in the inset in Fig. Elja), where we have 
assumed the following empirical^® magnetic field depen- 
dence of the gap: A(B)/A(0) = (1 - (B/BcY-^y-^. For 
the given parameters of our sample, Vt^e (0) = 2.4 mV and 
Vt,2e — 1.2 mV. Background charges will, however, mod- 
ify these thresholds. For an array of this size, Vt.2e and 
the part of Vt^e that does not depend on A become ap- 
proximately three times larger. At low field, we there- 
fore expect Cooper pairs to dominate, which is what we 
do observe for / > 100 fA. At larger fields, A is sup- 
pressed, and therefore also Vt,e, why at a given field the 
single-electron transport becomes significant. 

Inside the array, there is also a possibility that a 
Cooper-pair soliton centered on one island decays into 



two single-electron solitons centered on adjacent islands. 
This happens when the difference between the Cooper- 
pair and single-electron charging energies is smaller than 
2 A. Again disregarding background charges, the condi- 
tion for this is 

A(B)<;-^[l-exp(-l/A)], (2) 

which is satisfied for > 400 mT, meaning that we 
should detect pure le transport for higher fields. This 
agrees qualitatively with the data in Fig. [Sj^a), however, 
we note that thermal fiuctuations can break a metastable 
Cooper pair at lower fields. 

The arguments of the preceding paragraphs explain 
qualitatively the magnetic field dependence of the le — 2e 
crossover at relatively large currents. Let us now turn 
to the current dependent crossover that occurs for / < 
200 f A and _B| | < 350 niT and consider the different elec- 
tron and Cooper-pair tunneling rates. At the temper- 
ature of our experiment there should be practically no 
thermally excited quasiparticles, but experiments have 
shown that there are often non-equilibrium quasiparti- 
cles residing in the leads The threshold voltage for 
them to enter into the array is lower than that for Cooper 
pairs, 14,eNE = ^i,e(A = 0) = I4,2e/2. This means that 
for 14,eNE < V < Vt^2e only the non-equilibrium quasi- 
particles will enter and we should see pure le transport. 
Above Vt^2e we expect to see mixed le and 2e transport, 
and the different tunneling rates compete. The quasi- 
particle tunneling rate is proportional to the number of 
non-equilibrium quasiparticles present in the leads. This 
number is, in turn, determined through the competition 
between the process generating the quasiparticles, their 
recombination, and the le current, draining the quasi- 
particles into the array. The situation is similar to that 
of quasiparticle poisoning in the Cooper pair box, where 
a similar phenomenon has been observed^l and theoreti- 
cally described. This picture qualitatively explains the 
behavior seen in Fig. ^h) . An interesting aspect of this 
observation is that it should be possible to extract infor- 
mation about the density of non equilibrium quasiparti- 
cles by making more elaborate experiments of this kind. 

In conclusion, we have demonstrated time correlated 
tunneling of both individual electrons and individual 
Cooper pairs, and coexistence of the two, in a ID ar- 
ray of small Josephson junctions. We have shown that 
there is a crossover from single-electron transport to sin- 
gle Cooper-pair transport as a function of both the ex- 
ternal magnetic field and the current through the array. 
We describe the transport in terms of different thresh- 
old voltages for injection of charge into the array, and 
instability of Cooper pairs inside the array. 

We made the sample in the MC2 Nanofabrication Lab- 
oratory at Chalmers. We acknowledge helpful discussions 
with A. Kack, A. Zorin, and the members of the Quan- 
tum Device Physics Laboratory at MC2. This work was 
supported by the Swedish SSF, VR and by the Wallen- 
berg Foundation. 
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